Time-reversal symmetry breaking at Josephson tunnel junctions of purely d-wave 

superconductors 
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We study spontaneous time-reversal symmetry breaking at Josephson tunnel junctions of d-wave 
superconductors in the absence of subdominant components of the order parameter. For tunnel 
junctions, when the orientation is close to 0/45 (for which a gap lobe points towards the junction 
on one side and a gap node on the other), the mechanism of the symmetry breaking is the splitting 
of midgap states (MGS) by spontaneous establishment of a phase difference — ±7r/2 across the 
junction. This occurs for transparencies D 3> ^o/A and temperatures ksT <^ DAq, where ^0 is the 
coherence length, A is the penetration depth, and Aq is the maximum energy gap. On the other 
hand, tunnel junctions with D <^ ^o/X effectively behave as surfaces, for which the mechanism of 
symmetry breaking is self-induced Doppler shifts of MGS. For this instability, we calculate the phase 
transition temperature kBTrRSB = (l/6)(^o/A)Ao and show that the spatial shape of the gap is 
unimportant. 

PACS numbers; 74.50. +r, 74.20.-z 



The possibility of spontaneous time-reversal symme- 
try breaking (TRSB) states at various surfaces and in- 
terfaces of d-wave superconductors-has been intensively 
studied during the last few years.EI Some experimental 
findings, such as fractional flux quanta at certain grain 
boundariesB'0 and splitting pifnthe zero-bias conductance 
peak in zero magnetic fieldJj'B have been interpreted as 
realizations of this state. Theoretically, several different 
systems have been under consideration' surfaces p'0 twin 
boundariesjj and Josephson junctions0il3 with special 
orientations. 

In this paper we will study TRSB at Josephson tunnel 
junctions with orientation close to = 0, an = 7r/4 
(Fig. [1]) for arbitrary transparency D of the tunnel bar- 
rier. We emphasize that we are considering purely d-wave 
superconductors, meaning that a sub-dominant compo- 
nent of the order parameter is assumed to be absent. As 
will be shown, the TRSB effect is due to the specific prop- 
erties of the rpidgap states formed in these structures. 

In 1994 HuO showed that surface states with zero en- 
ergy, so-called midgap states (MGS), are formed at sur- 
faces and interfaces of c?-wave superconductors if the ori- 
entation angle a is non-zero. The largest spectral weight 
of the MGS appears when a d-wave gap node points di- 
rectly towards the surface/interface {a = 7r/4). 

It has been pointed out (see e.g. Refs. p],p^ that the 
large density of states exactly at the Fermi level associ- 
ated with the MGS is energetically unfavorable: if there 
exist mechanisms able to shift the MGS and produce a 
gap in the spectrum, the energy will be lowered and a 
phase transition into a state with broken time-reversal 
symmetry will take place. Splitting of MGS due to a 
complex d + is order parameter, with a subdominant 
surface/interface s-wave component, has been considered 



both at free surfacesQu and at Josephson junctionsO'E£l 
The possibility of instabilities at Josephson junctions of 
purely d-wave superconductors was first pointed out by 
Yipcl for the weak link case (transparency D = 1, no 
backscattering in the junction, see also Refs. [ll|Jl^). The 



TRSB state was later shown to be favorable also fpc 
nite but rather high transmissivity of the junction,t3"t3 
0.3 ^ D < 1. According to the symmetry arguments 
presented in Refs. p|,p"3[~p5|, the effect of TRSB in pure 
d-wave junctions heavily relies upon the non-sinusoidal 
current-phase relation in transparent weak links. Since 
the higher harmonics in the current-phase relation disap- 
pear in the tunnel limit, TRSB in tunnel junctions of pure 
d-wave superconductors is also expected to disappear 
and it was supposed that TRSB can only occur under 
such circumstances when a complex order parameter is 
formed.EJll3 However, it turns out that the symmetry ar- 
gument holds only for continuum states and for finite en- 
ergy bound states, but not for the MGS which contribute 
to the Josephson current with a term jpmportional to 
the first power in D at low temperature.E30 In the low- 
transparency limit it is not the non-sinusoidal current- 
phase relation itself that is important, but rather the un- 
even occupation of split MGS. Here we will show that this 
leads to TRSB also in tunnel junctions of purely d-wave 
superconductors. The driving mechanism of the instabil- 
ity is the displacement of MGS induced by spontaneous 
establishment of a finite phase difference 4> = ±7r/2 across 
the junctkm. ^similar to what happens in transparent 
j unctions Jj^E^"E3 In the extreme tunnel limit, the mecha- 
nism crosses over to self-sustained Doppler shifts q£MGS, 
which also produces instabilities at free surfaces .EfrEjEd 
Although the MGS contribution does not appear in 
conventional tunnel model calculations (despite the fact 
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that it is proportional to the first power in transparency 
D), the quasiclassical Green's function technique, in prin- 
ciple, includes it. However, since the necessary condition 
for TRSB for low transparency is ksT ^ DAq (corre- 
sponding to uneven occupation of split MGS, see below), 
the specific window of transparencies favoring TRSB 
found in Refs. [l^ , |l5| , 0.3 ^ D < 1, was actually a result 
of the choice of temperature, fc^T = 0.2kBTc ~ O.IAq. 

Consider now the Josephson tunnel junction in Fig. |^. 
We model the junction between the two clean two- 
dimensional d-wave superconductors by a square specu- 
lar barrier. In this case, the quasiparticle wavefunctions 
can be labeled by the conserved wave vector component 
parallel to the surface, ky = kpsinO, where kp is the 
Fermi wave vector. The gap functions in the supercon- 
ductors are 1\l/r{.x,9) = AogL/Fi{.x)cos[2{9 - ul/r)], 
where all angles are measured relative to the surface nor- 
mal (Fig. The fact that the gap may be suppressed 
near surfaces and interfaces of d-wave superconductors is 
reflected in the x-dependent functions gL/n^x). 

For the calculation of the dc Josephson current, one 
needs to consider contributions both from continuum 
states and Andreev bound states. If the two supercon- 
ductors were decoupled (zero transparency), there would 
be midgap surface states on the right side (a^ — n/A) 
and also finite energy surface states since the gap is sup- 
pressed near the surface forming a quantum well. In ad- 
dition there are gap edge states on the left side (a^ =0). 
For finite transparency the surface states form states of 
the entire junction. The energy of these Andreev states 
are shifted relative to the surface levels, the shift de- 
pending on the transparency of the barrier and the phase 
difference 4> across the junction. To clearly see the mech- 
anism of the TRSB instability we first consider a step- 
function dependence of the gap, gL/nix) = B(=Fa;). By 
solving the quasiclassical Bogoliubov-de Gennes equation 
for the junction, we find the energy of the midgap state 
(the ^-dependence of the order parameter is not explicitly 
written out here) 
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as plotted in Fig. ||(a). When a phase difference is ap- 
plied across the junction, the degeneracy of the ±fcj, 
MGS is lifted, as emphasized by the solid and dashed 
lines in Fig. ^(a). The contribution to the dc Joseph- 
son current from Andreev states is found via the relation 
Ix = {2e/fi){dE / d4))np{E) , where np the Fermi distri- 
bution function. The current carried by the bound states 
in Eq. (^ is therefore 
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where T is the temperature, Rn = Trh/e^kpLyD is the 
normal state resistance of the junction, Ly is the junction 
width, 7j = sm0, and D — J dr]D{r])/2. From Eq. (||) we 
find that the MGS contribution to the Josephson current 
at low temperature, /csT ^ |_Bo| ^ DAq, is of first order 
in the transparency D. In what concerns the continuum, 
zizky states are degenerate and carry current in opposite 
directions, which results in a cancellation of the main (of 
order D) current. The residual continuum contribution 
is small, of order D^. This cancellation is due to the 
sign change of A^; when we let ky ~ky, see Fig. |^. 
Both ky and —ky contributions have sin0 dependences, 
but they are 7r-shifted relative to each other (the sign of 
Ar is equivalent to a phase tt) and carry current ui-op^ 
posite directions. This is the symmetry argumenttHiJ'cj 
referred to above. Also the ztky bound states near the 
gap edges carry current in opposite directions. Although 
they are not degenerate (like the MGS they split under 
phase bias) they are equally populated for kgT <^ Aq 
and the sum of ztky currents is of order D^. We would 
like to emphasize that the MGS do not obey the ztky 
symmetry at low temperature because MGS with oppo- 
site signs of ky disperse on opposite sides of the Fermi 
level and are unequally populated. In Fig||(b) we plot 
the total current including the dominant MGS contri- 
bution and the small contributions from the gap edge 
states and the c.nijitiT][ii].m..sl.ates. In agreement with pre- 
vious workj3nl3liiiii3'E3aO the current-phase relation is 
TT-periodic. For increasing temperature, the MGS contri- 
bution is reduced, and when ksT ^ -DAq [the T = O.lTc 
curve in Fig. ^b)] the current is small, of order D^. In 
the intermediate region, the current in Eq. (||) scales as 
/Trlicn The above arguments leading to the domi- 
nating, of order D, MGS currents and small, of order 
D"^, non-MGS currents hold also for general spatial de- 
pendences of the gap functions. However, as shown in 
Ref. 23 , the numerical prefactor of the current calculated 
for step- function gaps is overestimated by a factor about 
two, due to an overestimation of the shift of the MGS 
with phase difference. However, this will not influence 
the instability we discuss in the following— 

By a phenomenological argument YipEI showed (see 
also the paper by Ostlundt3) that when the ±ky sym- 
metry cancels the first harmonic of the current in purely 
d-wave junctions with orientation ap ~ 0, = tt/A, 
time-reversal symmetry is broken if the parameters of 
the theory are chosen in such a way that the coefficient 
in front of the second harmonic is negative. The cur- 
rent due to MGS has indeed this negative sign. As a 
consequence, the equilibrium phase difference across the 
junction is 4'eq = ±7r/2, since the Josephson energy min- 
imum appears where the current through the junction is 
zero and the slope of the current-phase relation is pos- 
itive. That these phase differences really correspond to 
Josephson energy minima can be understood by noting 
that the energy of the midgap state is the lowest for 
(j)eq = ±7r/2, see Fig. ||(a) and also Eq. (||). Gonsidering 
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low temperature, when only the negative energy states 
are occupied, we see that the degenerate (peg = ±7r/2 
junction states correspond to occupation of ±ky time- 
reversed MGS, see the illustration in Fig. Assuming 
that the system chooses one minimum (0eg — 7r/2) or 
the other (— 7r/2), surface currents in the positive or neg- 
ative y-direction will appear and time-reversal symmetry 
is broken. The spontaneous surface current, calculated 
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where are the wavefunctions satisfying the quasiclas- 
sical Bogoliubov-deGennes equation, and np is the Fermi 
distribution function, is dominated by the MGS for the 
same reasons as the MGS dominate the Josephson cur- 
rent. In Fig. H(c) we plot jy{(f>) calculated to the right 
of the barrier at a; = 0. The surface current approaches 
its maximum value at the equilibrium phase differences 
(/>eg = ±7r/2, but it has opposite signs since MGS with 
opposite signs of ky are occupied at these two phase dif- 
ferences. 

Since the MGS surface current produces a magnetic 
field which costs energy, we must include this effect into 
the discussion of the instability. The spatial dependence 
of the magnetic field h is determined by the counter- 
flowing screening currents, and can be calculated via 
the superfluid momentum p^, h = — (c/e)V x p^. For 
type II superconductors, like the high- Tc superconduc- 
tors, Ps is found via the London equation V^Ps—Ps/A^ = 
{4:T:e/c^)jf^^{x)y = f{x)y. Because the MGS surface 
current remains finite in the limit D = 0, we are allowed 
to let 13 ^ and perform the calculation of the source 
current with the free surface MGS wavefunction. This 
allows us to explicitly take into account the spatial de- 
pendence of the gap. The surface current then takes the 
form 
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where C = J dl\A{l)\/hvF, x = - Ty^, and N = 
dxe~^'' is the normalization constant of the MGS 
wavefunction. The integration over trajectory angles rj 
is effectively cut off at 7y < 1 because of the tunneling 
cone described by D{rf). Since the MGS source only has 
a y-component, the non-trivial component of p^ is the y- 
component which only depends on x, Ps = ps{x)y. The 
solution of the differential equation, satisfying the bound- 
ary condition Unix— ,00 hz{x) = is Ps{x) — boe~^^^ + 
A c?a;'/(a;') sinh [(a: — a;')/A]. The constant 60 is fixed 
by the boundary condition hz{x = 0) = at the junc- 
tion, 60 = —A dxf{x)cosh{x/X). 

An important aspect of the screening problem is the 
separation of length scales: the surface current due to the 
MGS flows within a thin layer of width ^0 near the sur- 
face, while the screening currents flow in a much thicker 



layer of width A. For high-Tc superconductors the ra- 
tio $o/A ^ 1, and all quantities may be expanded in 
this small parameter. For this reason the convergence 
of the integrals in the expression for p^ is governed by 
the function f{x') which decays on the <^o length scale. 
Expanding the hyperbolic functions, the spatial shape of 
the gap function is cancelled in the leading term of the 
expression for ps : it appears both in the normalization N 
of the MGS wavefunction and in the integrals over f{x') 
which are integrals over the MGS wavefunction. Thus, 
the detailed spatial shape of the gap drops out of the 
calculation and the final form of ps is 
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We are now able to quantitatively study the difference 
in the thermodynamic potential f2 of junctions with and 
without broken symmetry: it consists of two parts, the 
energy cost of having a magnetic field and the energy 
gain due to the shifts of MGS 
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where the second line is valid in the low-temperature 
limit keT < l^ol ~ -DAn. Clearly, for D > ^o/A, 
Ari < and there is an instability. 

Rotating the superconductors away from the = 0, 
an = 7r/4 orientation, the equilibrium phase difference 
across the junction is shifted continuously away from 
±7r/2 towards either or ±7r depending on the direction 
of rotation. This happens because the ±fcy symmetry is 
lost and non-MGS contributions are able to dominate. 
However, a numerical calculation shows that in the low- 
temperature region, ksT <C -DAq, where the MGS dom- 
inate, the TRSB, that is 7^ or ±7r, is quite robust: 
it survives rotations up to 10°. 

From Eq. ^ it is found that when D < ^o/A the 
energy of the magnetic field may become larger than 
the Josephson energy and the instability is lost. How- 
ever, for such small transparencies it is necessary to take 
into account the Doppler shift of the MGS due to the 
finite superfluid momentum, which will assist and up- 
hold the TRSB instability. In the low transparency limit 
D -C $o/A, for which the Doppler shift is much larger 
than the shift due to finite phase difference, the junction 
behaves like a free surface with the MGS localized on the 
right side of the junction. Using the above calculation but 
taking the energy of the MGS to be Eq{?]) = ■ vp — 
PsVpT], one finds from Eq. (^) in the zero-temperature 
limit An = [Ep /16\){1 - 4)[1 + 0(Co/A)] < show- 
ing that the TRSB state is favorable. Eq. (||) with 
EqW) = PsVpTj takes the form of a self-consistency equa- 
tion for Ps(T). Near the phase transition temperature. 
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Ps is small and the inequality PsVp << ^ksT is ful- 
filled; this allows an expansion of the hyperbolic function, 
which leads to \ps\ — 2yy5/3{h/ X)^^! — T/Ttb.sb, where 
ksTTRSB = (1/6) (Co/A) Ao is the temperature where the 
second order phase transition into the TRSB state occurs^ 
This surface instability was discovered by HigashitaniEEl 
in connection with a study of the paramagnetic response 
of the MGS to external magnetic fields. Our present 
calculation shows that the surface instability is not sen- 
sitive to the spatial profile of the gap; since we consid- 
ered a finite value of the parameter ^o/ A we rwere able 
to calculate the phase transition temperaturccJ Within 
the framework of the quasiclassical Green's function tech- 
nique-analogous results were recently obtained by Barash 
et alEB 

For orientations of the surface different from a = 7r/4, 
MGS exists only for trajectories satisfying sgn(AA) = 
— 1. The surface current is then reduced, leading to 
a continuous reduction of the phase transition temper- 
ature, ksTTRSB = (l/6)(Co/A)Ao (cos^fo-sin^fo), 
with misorientation 5a = |7r/4 — a| € [0,7r/4]. For small 
misorientations, <C 1, the reduction is quadratic in 
5a and the instability is not dramatically sensitive to the 
exact orientation of the surface. Note, however, tiia^he 
instability is quite sensitive to surface roughness .E§E3 

In conclusion, we have studied time-reversal symme- 
try breaking at Josephson tunnel junctions of d-wave 
superconductors assuming that a subdominant compo- 
nent of the order parameter is absent. For the D ^ 
Co/A junction case, at temperatures ksT <^ DI^q for 
which MGS contribute to the Josephson current with a 
term proportional to the first power in _D, the TRSB is 
due to spontaneous establishment of a phase difference 
(f) = ±7r/2 across the junction which splits the MGS and 
produces a surface current. In the extreme low trans- 
parency limit, -D <C Co/A, the mechanism responsible for 
the instability is instead the same as for the free sur- 
face. In this case, below a phase transition temperature 
ksTTRSB — (1/6) (Co/A) Aq, it is energetically favorable 
to have Doppler shifted MGS carrying a surface current 
and associated screening currents upholding the Doppler 
shifts. The detailed spatial shape of the gap does not 
influence this instability. 
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FIG. 1. The specular d-wavc Josephson junction under 
consideration: the right superconductor has a gap node point- 
ing towards the junction {aa = 7I-/4), while the left supercon- 
ductor has a gap lobe pointing towards the junction (qz, =0). 
For the surface problem we have vacuum for a; < 0. Shown 
are also the scattering events, consecutive normal scattering 
at the junction and Andreev reflection in the superconductor, 
leading to the formation of midgap states. 
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FIG. 2. (a) Dispersion of the Andreev bound state with the 
phase difference across the junction. Solid and dashed lines 

are bound states at angles 6 = n/9 and 6 = —n/9 respec- 
tively, D{9 — ±7r/9) ~ 0.01. (b) The current-phase relation 
at three difl'crcnt temperatures: T = O.lTc, T — O.OlTc, and 
T = O.OOlTc, for D w 0.009 {hrd = 5, U = 1.2Ef, where 
U and d are the height and width of the barrier and Ef is 
the Fermi energy). For decreasing temperature, the current 
crosses over from being of order to D and becomes in- 
creasingly non-sinusoidal, (c) Phase dependence of the sur- 
face current density calculated to the right of the barrier. 
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